PAMELA and dark matter 
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Assuming that the positron excess in PAMELA satellite data is a consequence of annihilations of 
cold dark matter, we consider from a model-independent perspective if the data show a preference 
for the spin of dark matter. We then perform a general analysis of annihilations into two-body 
states to determine what weighted combination of channels best describes the data. 



Introduction. One of the unsolved problems in cos- 
mology and particle physics is the nature of dark matter 
(DM) which accounts for about 20% of the energy den- 
sity of the universe. Particle physics models typically 
relate discrete symmetries with the existence of a sta- 
ble cold DM candidate. A variety of such models have 
been suggested that provide viable explanations of the 
DM. Weakly interacting massive particles (WIMPs) can 
be either of integer or noninteger spin. The classic case 
of supersymmctry (SUSY) has a spin- 1/2 neutralino as 
dark matter whereas extra dimensional models and col- 
lective symmetry breaking models have spin-1 dark mat- 
ter. Specific realizations are the minimal Universal Extra 
Dimensions (mUED) [lj and Little Higgs with T-parity 
(LHT) 2] models. Spin-0 dark matter is possible in mod- 
els with an additional singlet in the scalar sector of the 
Standard Model 0. 

Recent evidence for a positron excess in the Payload for 
Matter Antimatter Exploration and Light-nuclei Astro- 
physics (PAMELA) data [4] spurs attention to WIMPs 
whose annihilations in the galactic halo can explain an 
excess over backgrounds d, @ • PAMELA data presented 
thus far show a turn-up in the energy spectrum at about 
10 GeV and a steady rise up to 100 GeV with no fall- 
off in that dataset. 1 The shape of the spectrum bears 
directly on the annihilation mechanism. Spin-1 DM an- 
nihilations directly into e + e~ produce a line spectrum, 
whereas spin- 1/2 Majorana DM will give a continuum 
spectrum from secondary decays of weak bosons, quarks 
and lcptons. 

Our goals are to study in a model-independent ap- 
proach if the PAMELA excess provides hints about the 
spin of the DM particle and what annihilation channels 
are favored by the data. We do not subscribe to any 
specific particle physics model, but comment on models 
where appropriate. We also do not require that the mea- 
sured relic abundance be reproduced since this is highly 
model-dependent. Moreover, the total energy density in 
dark matter may be comprised of several components, 



We do not consider the excess in the e+ + e energy spectrum 
between 500 and 800 GeV seen by the PPB-BETS balloon ex- 
periment 0. 



so only an upper bound need be imposed on the energy 
density of a particular DM particle. The nature of our 
analysis precludes us from making projections for signa- 
tures at IceCube, the Fermi Gamma-ray Space Telescope, 
direct detection experiments and the Large Hadron Col- 
lider, all of which are interesting in their own right. 

Modelling the positron signal and background. 

The positron background expected primarily from super- 
novae and from collisions of cosmic ray protons and nuclei 
on the interstellar medium is simulated in Ref. [8]. The 
results of the simulation have the convenient parameteri- 
zation @, <S> b e k + 9 = 4.5£ a7 /(l + 650£ 2 ' 3 + 1500£ 4 - 2 ), with 
the energy of the positron E in GeV. Since we present 
our results as a positron fraction <E> e + / (<& e + + $ e - ) which 
allows for cancellations of systematic uncertainties and 
the effects of solar activity, we also need the electron 
background which is analogously parameterized by [To| 
^J> = Q.l&E- ll /{\ + HE 09 + 3.2E 215 ) + 0.7£°'7(1 + 
110-E 1 - 5 + 580-E 4 ' 2 ). Solar modulations arise from the 
phase of the solar cycle and from the opposite charges of 
electrons and positrons. Without charge sign bias, the 
positron ratio is independent of solar activity. However, 
since PAMELA data show evidence of charge sign de- 
pendence for positron energies below 10 GeV, we only 
analyze data above 10 GeV. 

Positrons produced in DM annihilations propagate 
through the interstellar medium to the earth and as a 
consequence suffer absorption effects that broaden the 
positron spectrum to lower energies. We estimate the 
primary positron flux from dark matter annihilations ac- 
cording to the prescription of Refs. [TH El, Here 
we briefly describe the procedure and refer the reader to 
Refs. [Ifl 01 for details. 

The positron number density per unit energy is gov- 
erned by the diffusion-loss equation with diffusion coef- 
ficient K(E) = K Q (E/GeV) s which describes propaga- 
tion through turbulent magnetic fields, and is taken to 
be independent of position. The equation also accounts 
for energy losses through synchrotron radiation and in- 
verse Compton scattering on the cosmic microwave back- 
ground and infrared galactic starlight. The diffusion zone 
in which the diffusion-loss equation is applicable is mod- 
elled as a cylinder of height 2L and radius 20 kpc that 
straddles the galactic plane in which most cosmic ray in- 
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Model 
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K (kpcVMyr) 


L (kpc) 


Min 


0.55 


0.00595 


1 


Med 


0.70 


0.0112 


4 


Max 


0.46 


0.0765 


15 



TABLE I: Three sets of parameters describing cosmic ray 
propagation [l^| . The Med set is the best- fit to the mea- 
sured boron to carbon ratio. The Min and Max sets minimize 
and maximize the positron fluxes above 10 GeV, respectively. 



teractions take place. The positron number density is 
assumed to vanish on the surface of the cylinder, since 
outside the diffusion zone the positrons propagate freely 
and escape into the intergalactic medium. The source of 
the positrons due to DM annhilations depends on the DM 
density profile and on the annihilation cross section. For 
the former, we adopt a cored isothermal halo profile [l3[. 

The normalization Kq and the spectral index S of the 
diffusion coefficient, and L can all be selected to be con- 
sistent with the measured boron to carbon ratio in cosmic 
rays [l4j. We consider three sets of these parameters, 
"Min" , "Med" and "Max" , of Ref. and reproduce 
them in Table |TJ The Med set has values of Kq , 5 and L 
which best fit the measured boron to carbon ratio. The 
Min and Max sets minimize and maximize the positron 
fluxes above about 10 GeV. Needless to say, the Min and 
Max sets are only representative, since the positron flux 
depends on the mass of the DM particle Mdm and on 
the annihilation channel once the halo profile is selected. 

Assuming steady-state conditions, a semi-analytic ex- 
pression for the primary positron flux has been ob- 
tained [ll|, [l2j • The result depends on a so-called "halo 
function" which encodes the physics of cosmic ray prop- 
agation. We employ numerical fit functions [T3| for the 
halo functions pertinent to the isothermal profile with 
the Min, Med and Max propagation parameter sets. We 
allow for the possibility of high density substructures in 
the dark matter halo that enhance the positron flux by 
an energy-independent "boost factor" B. Note that N- 
body simulations suggest that B can not be larger than 
about 10 and may be energy-dependent [151 ]. 

Dark matter annihilations and spin. To begin 
with, we assume Mdm is smaller than the top quark 
mass. This choice is dictated by our interest in model- 
independence. (We shall see later, by extending the 
range of Mom, that the positron data typically select 
DM lighter than the top quark for the Med set). If the tt 
channel were open, the relative contributions of different 
channels to the positron spectrum would depend on the 
details of the Higgs sector. 

Since the PAMELA data show a sharp rise, we ini- 
tially only consider positrons from a e + line spectrum or 
from the two-body decays of pair-produced weak bosons 
at the source. Specifically we study the spectra from 
direct production, DMDM — > e + e~ (which produces a 
positron line close to Mdm), and from secondary produc- 



spin 


s-channel 


t, it- channel 


t, M-channel 




Higgs 


fermion 


boson 





LL, TT 


X 


LL 


i 

2 





TT 


X 


1 


LL, TT 


X 


LL, TT 



TABLE II: Polarizations of W pairs produced by static an- 
nihilations DMDM — » 1P + W~ depend on the spin of the 
DM particle. "LL" and "TT" indicate that the W bosons are 
longitudinally and transversely polarized, respectively. "X" 
indicates that there is no contribution at the tree-level, and 
"0" indicates that the amplitude vanishes in the static limit. 
Note that Dirac fermion DM also has contributions from s- 
channel Z-exchange. 



tion from the process DMDM — > W + W~ ? Concrete 
examples of direct annihilation into e + e~ are found in 
mUED and LHT in which spin-1 DM annhilate by ex- 
change of an odd-parity fermion [ltil \v\ . Direct annihi- 
lation also occurs for hidden/mirror Dirac fermions and 
sterile neutrinos. The latter constitute warm DM which 
is not relevant to our study of nonrelativistic DM. If DM 
is a Majorana fermion, helicity suppression prevents the 
direct production of e + e~. For scalar DM the amplitude 
for static annihilation into light fermions vanishes [Lsl ]. 
Since the production of W pairs is spin-dependent, we 
further classify the positron spectra according to whether 
the W bosons are longitudinally polarized or tranversely 
polarized. 

The normalized distributions for the e + energy are 
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for the WW transverse (TT) modes and longitudinal 



(LL) mode, respectively, where 0^ = 1 



' w 



/Af 



DM 



and x = 2E e + /M DM <$1 ± Pw- In general, if the W + W~ 
channel has both TT and LL contributions of rela- 
tive weights a and b, then the resultant distribution 
which combines two TT modes and one LL mode is 
(af LL (x) + bf TT (x))/{a + b). 

In Table UH we categorize the polarization modes of 
the W pair according to the spin of the DM particle. 
While fcrmionic DM can not annihilate into W pairs 
via s-channel Higgs exchange in the static limit, spin- 
and spin-1 DM annhilations (with relative weights 
a = (1 + /3^/) 2 and 6 = 2(1- flw) 2 ) give the distribution, 



1 dN _ 3[1 + / 
N~dx ~ 2f3 w (3 
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2 Although each Z in a Z pair produces a positron, since 
a(DMDM -» W+W-)/cr(DMDM -> ZZ) rj 2 in the high- 
energy limit, and since the leptonic branching fraction for Ws 
is three times as much as for Zs, W pairs produce three times as 
many positrons as Z pairs with almost identical distributions. 
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FIG. 1: Annhilations of DM directly into e + e~ give the e + 
line at about Mom ~ 150 GeV. The secondary positron spec- 
trum from decays into e + v of longitudinal (transverse) W 
bosons is labeled WlWl (WtWt); the soft component of the 
spectra are neglected for illustration. Including the soft com- 
ponent (with spin-correlations neglected) results in a much 
softer spectrum labelled WW that does not fit the PAMELA 
data above 10 GeV. The solid curve is the expected back- 
ground. The Med set of propagation parameters is used with 
a cored isothermal profile for the DM halo. 
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FIG. 2: Positron fraction from DM annihilation into e + e~ , 
fi + , t + t~ , W + W~ , ZZ, cc, bb, and hh, with a Standard 
Model Higgs boson h of mass 120 GeV for the Med set of 
propagation parameters. We have assumed that the DM an- 
nihilates into each mode with a 100% branching fraction. Ac- 
counting for the smaller boost factor, the e + e~ mode is some- 
what preferred; see Table Hill The Max set yields spectra very 
similar to the Med set. 



Note that as [3w ~ * 1> the longitudinal mode dominates. 
The DM particle in both mUED and LHT can annihi- 
late via s-channcl Higgs exchange. Whether these mod- 
els produce line or continuum spectra or both depends 
on specific realizations. 

Fermionic DM annhilations via t- or u- exchange of a 
fcrmion give only TT modes. The positron spectrum is 
then simply /tt(^) [HI- SUSY provides the common 
example of neutralinos that annihilate dominantly by t- 
and u-channel chargino exchange. 

Analysis. In Fig. [TJ we show spectral distributions of 
positrons produced in annhilations of a DM particle of 
mass 150 GeV that fit the PAMELA excess. We have as- 
sumed that when direct production occurs, annihilations 
into W + W~ are negligibly small. For the W + W~ chan- 
nels, only the hard spectra from W + — > e + v are shown. 
The soft component of the spectra from the W other de- 
cay modes is neglected to emphasize the small difference 
between the hard spectra from WlWl and WtWt- In 
what follows, we disregard the effects of spin-correlations. 
Including the soft component, we find annihilations dom- 
inantly into W + W~ do not provide a satisfactory spec- 
trum. 

We now enlarge the scope of our study by allow- 
ing Mom to be as large as 1 TeV and allowing arbi- 
trary weights for the following annihilation modes: e + e~, 
t + t~, W + W~ , ZZ, cc, bb, ti and hh, with a 
Higgs boson h of mass 120 GeV which will decay primar- 
ily into b and r. Annihilations into Zh can be accounted 
for by the average of the ZZ and hh channels. The sub- 
sequent decays were computed using micrOMEGAs [2(| ■ 
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FIG. 3: Same as Fig. [2] but for the Min set of propagation pa- 
rameters. The /i + fi~ , t + t~ and W + W~ modes are preferred; 
see Table iHll 



We denote the weights by f xy , where for example, f e + e - 
is the weight of the e + e~ channel. In Figs. [2] and [31 we 
show the positron fraction from each of these channels 
(except tt) with f xy = 1 for Mom — 150 GeV. It is 
evident that annihilations into bosons and quarks yield 
too soft a spectrum, while annihilations into leptons are 
easily compatible with the data. The corresponding \ 2 
values are listed in Table IIIII We have not displayed re- 
sults for the Max set because it gives spectra similar to 
those for the Med set. Although the lowest x 2 P er de- 
gree of freedom is 2 for the Min set, we do not reject this 
parameter set given that uncertainties in the positron 
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TABLE III: \ 2 for positron spectra from two-body annihila- 
tions of DM with mass 150 GeV for the Med and Min models 
of cosmic ray propagation. Results for the Max set of param- 
eters are similar to those for the Med set. The number of 
degrees of freedom in each case is 6. 
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background have not yet been estimated. Since the Med 
set has a larger diffusion zone height 2L than the Min 
set, the flux of positrons incident at PAMELA is larger, 
thus requiring a smaller B. This explains the mode-by- 
mode lower boost factors for the Med set in Table IIIII 
On the other hand, the Min set has a smaller spectral 
index 5 with relatively weaker diffusion at higher energy, 
resulting in a harder positron spectrum. This is evident 
from a comparison of the spectra in Figs. [5] and [3J 

We perform a Markov Chain Monte Carlo (MCMC) 
analysis by varying Mjjm between 100 GeV and 1 TeV, 
the boost factor B and the weights f xy between and 1 
to determine the combination of annihilation modes that 
fits the positron data best; see Ref. [IH for a descrip- 
tion of our MCMC methodology. We set the annihilation 
cross section to be (av) = 3 x 10~ 26 cm~ 3 s _1 , which is 
the typical value required to reproduce the observed relic 
abundance barring co-annihilations with other particles. 
The probability distribution of / is shown in Fig. Q] For 
the Med set of propagation parameters, the e + e~ mode is 
preferred. For the Min set, a preference for the fi + fj,~~ and 
t + t~ modes is evident, and the e + e~ mode is not favored 
over the non-lepton modes. In general, the data show a 
preference for lepton modes. The correlation matrix for 
the nine modes is shown graphically in Fig. O The 2a 
contours in planes of weights for pairs of modes are plot- 
ted after marginalizing over all other modes. There is 
essentially no correlation between modes. It is notewor- 
thy that while the 2a region for the [i + fx" mode and any 
mode (other than e + e~) is consistent with (0,0) for the 
Med set, it is not so for the Min set. This is because a soft 
component is necessary to fit the data for the Min set. 
The probability distribution of B is shown in Fig. [6j For 
the WIMP annihilation cross section we have adopted, 
the boost factor is about 50 for the Med set, which is not 
grossly unreasonable. 

From Fig. [Jj we see that the range of DM masses fa- 
vored by the positron data depends on the details of cos- 
mic ray propagation. At 2a, Mdm is below 215 GeV for 
the Med set and below 445 GeV for the Min set. The 



FIG. 4: The probability distribution of the weight of each 
channel / for the Med and Min sets. The medians, and la 
and 2a C. L. ranges are indicated in the lower panels. 
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FIG. 5: A graphical representation of the correlations between 
modes for the Med (upper triangle) and Min (lower triangle) 
propagation sets. The cells along the diagonal show the prob- 
ability distribution of / corresponding to the mode labeled. 
The contour plots show the 2a allowed regions in planes of 
weights, with / of the column (row) mode along the x-axis 
(y-axis). No two modes are significantly correlated with each 
other. 



correlation between B and Mdm in Fig. [8] shows that 
ligher DM particles require a smaller boost factor to ex- 
plain the PAMELA positron excess. Also, with the Med 
set of propagation parameters, very large boost factors 
are avoided. 

The antiproton spectrum measured by PAMELA up 
to 100 GeV 1211] shows no deviation from the expected 
background [22j (which has larger uncertainties than the 
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FIG. 6: The probability distribution of the boost factor B for 
the Med and Min propagation sets. 




100 150 200 250 300 350 400 200 400 600 800 1000 

M DM M DM 



FIG. 7: The probability distribution of M DM for the Med 
and Min sets. The medians which are 125 GeV (Med) and 
170 GeV (Min), and la and 2a C. L. ranges are indicated in 
the lower panels. 



positron background associated in part with the consider- 
ably greater propagation distance of antiprotons). Since 
our approach is model-independent we can not make def- 
inite statements about consistency with the cosmic an- 
tiproton data. By choosing an appropriate boost factor 
(which can be different from that for positrons) and ap- 
propriately modelling the propagation of antiprotons, it 
is easy to remain in agreement with the data. Within 
our approach it is also possible to have consistency by 
suppressing the DM annihilation branching fraction to 
antiquarks. As an illustration, in Fig. [9] we show the 
antiproton to proton flux ratio measured by PAMELA, 
and the theoretical expectation for the W + W~ channel 
from annihilations of DM of mass 150 GeV. Boost fac- 
tors for the antiproton flux below 3.3 yield agreement at 
the 2a C. L. The light dashed curve shows the p/p flux 
ratio if the antiproton boost factor is taken to be equal 
to the positron boost factor that fits the positron spec- 
trum. Clearly, different boost factors are necessary. The 
1-2 orders of magnitude difference in the e + and p boost 
factors is a problem. 



FIG. 8: The correlation between B and Mom for the Med 
and Min propagation sets. 
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FIG. 9: The p/p flux ratio measured by PAMELA is consis- 
tent with the expected background (solid). The dark dashed 
curve is the expected spectrum from DM annihilations to 
W + W~ , allowing for a boost factor (equal to 1) that is dif- 
ferent from that for positrons. The light dashed curve shows 
the spectrum if the boost factor for the positron fraction is 
applied to antiprotons; the inset shows a magnified view. 



Conclusions. Our results are summarized in the fig- 
ures. We have shown that the PAMELA positron excess 
does not favor a DM particle of a particular spin. The 
data do not discriminate between positron spectra from 
direct production and from secondary decays of polar- 
ized W bosons. However, PAMELA is expected to col- 
lect positrons up to about 270 GeV. With those data 
it should be possible to draw stronger conclusions. If 
the data show a line, popular SUSY models will be in 
danger of being excluded and models with extra dimen- 
sions and collective symmetry breaking will gain support 
since they have spin-1 DM. Models with Dirac fermions 
as DM will also be viable. On the other hand if the data 
roll-over smoothly near the endpoint, and are fit well by 
positrons from transversely polarized W bosons, SUSY 
will be indicated. If positrons from longitudinally polar- 
ized W bosons are preferred by the data, neutralino DM 
will be in jeopardy, and the DM candidates of mUED 
and LHT will be preferred. To make such fine distinc- 
tions in spectral shapes will require much larger datasets 



6 



from PAMELA and the Alpha Magnetic Spectrometer. 

By considering nine different two-body annihilation 
channels with arbitrary weights, for dark matter lighter 
than 1 TeV, we found that lepton modes are generally 
preferred by the positron data, and which lepton modes 
are favored depends on the details of cosmic ray propa- 
gation. The /i + [i~ and t + t~ modes fit the data better 
than the obvious e + e~ mode for the Min set. Also, we 
found that dark matter masses selected by the data de- 
pend on the propagation model. The 2a upper limit is 
215 GeV for the Med set of propagation parameters and 
445 GeV for the Min set. Results for the Max set are 
similar to those for the Mid set. 



It is important to bear in mind that although astro- 
physical processes are expected to produce a positron 
background that falls with energy, it may still be that 
astrophysical sources such as pulsars could mimic the pu- 
tative DM signal. Confidence in the DM interpretation 
will be strengthened by signals in other experiments, in- 
volving both direct and indirect detection methods. 
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